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Abstract 

We discuss the space-time determinantal random field which arises 
for the PNG model in one dimension and resembles the one for Dyson's 
Brownian motion. The information of interest for growth processes is 
carried by the edge statistics of the random field and therefore their 
universal scaling is related to the edge properties of Gaussian multi- 
matrix models. 

1 Introduction 

There is a huge variety of stochastic growth processes and we refer to the 
recent monographs by Barabasi, Stanley[l] and Meakin[2]. Amongst them 
the KPZ (Kardar, Parisi, Zhang) growth models enjoy a particular popularity 
in theoretical circles. Prototypes are the Eden growth, where each perimeter 
site of the current cluster is filled after an exponentially distributed waiting 
time, and ballistic deposition, where there is a low intensity random flux of 
incoming particles which then attach to the current surface profile. Thus in 
broad terms KPZ growth is characterised by being stochastic with a 

• local growth rule, 

• smoothening mechanism. 
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The list of experiments well described in terms of KPZ is rather short, see 
Mylls et al. 3 for recent experiments on slow combustion fronts. But one 
obvious reason for its popularity is that KPZ growth is rather close to the 
stochastic dynamics commonly studied in Statistical Mechanics with the fine 
twist that it does not satisfy detailed balance (KPZ growth is not stochas- 
tically reversible). There is a second reason, however: through a Cole-Hopf 
type transformation KPZ growth maps to directed first passage percolation. 
Thereby techniques and insights from the theory of disordered systems come 
into play. In fact, there are some rather close analogies, one example being 
the rigorous discussion of the overlap for directed polymers [I]. 

Just before the previous ICMP congress it came as a total surprise that 
KPZ growth in one spatial dimension (the surface is the graph of a function 
over R) is linked to Gaussian random matrices. For example, the height, 
h(r), at time r above a given reference point grows at constant speed t> 
with some random fluctuations, 

h(T)=V T + T X '*Z (1) 

for large r. If at the reference point the macroscopic profile has a non-zero 
curvature, then £ has the same distribution function as the largest eigen- 
value of a GUE random matrix, which in the community is known as the 
Tracy- Widom distribution function F 2 , P(£ < x) — F 2 (x). F 2 is related to 
the Hastings-McLeod solution of the Painleve II differential equation. The 
scaling exponent 1/3 and the random amplitude ^ in (|T|) are expected to be 
valid for all one-dimensional KPZ growth models. 

The purpose of our contribution is to explain how this connection arises. 
The rough indication can be given already now: For some very particular 
growth models in the KPZ class, there is a determinantal process in the 
background. Its structure has some similarity to the determinantal process 
appearing in the context of Gaussian multi-matrix models. Since the relevant 
information is linked to edge scaling, in the scaling limit the GUE edge 
distribution arises. 

2 GUE, Dyson's Brownian motion, and the 
Airy process 

Dyson [3] considers an Ornstein-Uhlenbeck process, A(t), t G JR., on the space 
of N x N Hermitian matrices. Its transition probability is given through a 
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Mehler type formula as 

F(A(t) E dA\A(0) = A) = ± exp[-iTr(A' - qAf/{\ - q 2 )]dA' (2) 

with q = e~*. In particular, for t — ► oo the process converges to the GUE 
ensemble given by 

Z-'expl-^TiA^dA. (3) 

For the joint distribution of the stationary process at times ordered as to < 
ti < ... < t m one obtains therefore 

F({A(t J )edA J ,j = 0,l,...,m}) 



z exp 



m m 

X ( E Tr (^' - ^Vi)V(l - <$) + Tr^) ] J] dA 3 (4) 

3=1 3=0 



with qj = exp[— (tj — is the Gaussian multi-matrix model of the 

title. 

The process A(t) induces a process on the eigenvalues Ajv(t) < ••• < Ai(£) 
of A(t). It is stationary by construction and happens to be again Markov. 
In fact, it satisfies the set of stochastic differential equations 
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d\ j {t) = (~\ j (t)+(p/2) Mtymy^dt+db^t), j = i,...,n, 

(5) 

with {bj(t), j = 1, N} a collection of N independent standard Brownian 
motions. In case of Hermitian random matrices (3 = 2. The repulsive drift 
ensures that A 3 -+i(t) < Xj(t) for all t, any j. Clearly (JHJ) has the unique 
stationary distribution 



z exp 



-. N N 

4e\ 2 1 n w-^u^- w 



iV 

j=l l<i<j<N j=l 



There is a way to rewrite (0) which turns out to be computationally very 
powerful. We define a random field <(>n over M? through 
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Mx,t) =^2S(x-X j (t)). (7) 



Then 0at is determinantal in the sense that its moments (correlation func- 
tions) have a determinantal structure of the form 

in 

(JJ <Pn(xj, tj)) = det{R N (xi, ti\ xj, tj)}™ j=1 (8) 

3=1 

for distinct times t\, . . . , t m . The defining kernel Rn can be written in terms 
of the harmonic oscillator Hamiltonian 

H » = -\% + ^fS*> ( 9 ) 

which has eigenvalues E n — n/N, n — 1, 2, ... . Let -fTjv be the Hermite kernel 
which is the spectral projection onto {Hn < 1}. Then 

R N (x,t;x',t') = (e- tHN (K N - lQ{t - t'))e t,HN yx,x') (10) 

with 0(t) = 1 for t > and B(t) = for £ < 0. In particular, at t — (or 
any other time by stationarity) 

m 

(J] cf>( Xj ,0)) = det{^(a; i , % )}-. =1 , (11) 

3=1 

as well known from the random matrix bible by Mehta[6j. 

A natural question is to study the statistics of lines close to the, say 
upper, edge. In our units the top line fluctuates at level \/2N, as can be seen 
by equating the Fermi energy E-p, E-p = 1, with the energy of the confining 
potential of ©. Thus we shift our attention to x = \/2N and linearise there 
the potential, a procedure which should be accurate for large N. Then the 
imaginary time Schrodinger equation for © goes over to 

d t ^ = {- l -d 2 x + ^V2x)^. (12) 

It becomes iV-independent under the scaling t N 2 ^ 3 t, x ~-> iV 1//3 x/v^2, 
resulting in the Schrodinger equation with Airy Hamiltonian 

d t ^j = H^, H = -d 2 x +x. (13) 
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We have identified the edge scaling and conclude that, in the sense of 
convergence of moments, 



lim -=N 1/3 (f> N (V2N + -^=N 1/3 x , N 2/ H) = <j>(x, t) . (14) 
n^oo V2 

The prefactor comes from the spatial volume element when integrating both 
sides in (fUjl over a compactly supported test function. is called the Airy 
random field. Since 4>n is determinantal, so must be its limit. Hence, for 
distinct times t 1 , ...,t m , 



(II ^> = det {^, U- Xj , ' ( 15 ) 

compare with (|SJ). The defining kernel is now given through 
R(x,t;x',lf) = (e- tH {K -m{t-t'))e t,H Yx,x') (16) 
= Bign(f - 1) [ dXQ(X(t - t'))e A( *'- 4) Ai(a; - A)Ai(x' - A) , 



where sign(t) = 1 for t > and sign(t) = — 1 for t < 0, Ai the Airy function, 
and K the spectral projection onto {H < 0}. K is known as Airy kernel. 
For fixed t the limit is studied by Forrester [7], and Tracy, Widom[8j. Some 
aspects of the multi-matrix model are discussed by EynardjH]. 

The top line of the Airy random field is called the Airy process ^U]> de- 
noted by A(t). Its joint distributions can be written in a concise way. Let 
t\ < ... < t m and let us define the operator, R^ m \ on L 2 (M x {1, ...,m},dx) 
through the integral kernel R(x, tf, x', tj)x({x > £i})x({ x> > hJ = 

1, ...,m. Then R^ is trace class and 

P(M(*i) <Zi,-,A{t m ) <U}) = P({#M;)=0 for x>£ j ,j = l,...,m}) 

= det(l - R {m) ) . (17) 

The Airy process has continuous sample paths and is stationary, by con- 
struction. Some more explicit expressions for joint distributions are given in 
[HI H2J HSJ HI]. In particular, (^l(O) 2 ) = 0.81325..., ((.4(0) - A{t)f) = 2\t\ 
for small t, and {A(0)A(t)) - (-A(O)) 2 = t~ 2 + C(t" 4 ) for large t. 
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3 Polynuclear growth, its determinantal ran- 
dom field, and edge scaling 

For the polynuclear growth (PNG) model the height function at time r takes 
integer values, h(x, r) G Z, x G R, r > 0. xh r) has jumps of unit size 
only. Upward steps move, in the x direction, with velocity —1 and downward 
steps with velocity 1. They annihilate each other upon collision, which is the 
required smoothening mechanism. The surface grows through nucleation. 
At such an event, say (xi,Ti), the height h(x, T\) is increased by one unit 
at Xi, thereby creating an upward and downward step, which move apart 
immediately under the deterministic part of the dynamics. The nucleation 
events are Poisson in space-time with intensity 2. We consider the droplet 
geometry which is enforced by allowing nucleations only in the space interval 
[ — t, t] . Initially h(x, 0) = and, by assumption, h(x,r) = for \x\ > r. 
Obviously, according to our criteria of the Introduction, the PNG model is 
in the KPZ class. 

For PNG, as just described, there is no determinantal process in sight. 
The miracle happens through the RSK (Robinson, Schensted, Knuth) con- 
struction. The idea is to extend the model with additional lines which record 
the information lost in annihilation events. Thus we introduce the lines 
hj(x,r), j = 0,-1,-2,..., and set h (x,r) = h(x,r). Initially hj(x, 0) = j. 
h evolves according to the PNG specified above. Given 1i (x,t) all lower 
lying lines have a deterministic dynamics. The steps of line hj, j < — 1, 
move with velocity ±1, as before. A nucleation, say at space-time point 
(x,t), takes place whenever in line j + 1 an upward and downward step an- 
nihilate each other at (x,t). Thus x h- > hj(x,r) has jumps of unit size only, 
hj(x,r) = j for \x\ > r, and hj(x,r) > hj-^x, r). Furthermore there is a 
random index j such that, for j < j , hj(x, r) = j for all x. Let us set, at 
fixed time r, 



with j G Z, t G [— r, r]. i] T (j,t) is a determinantal random field over Z x 
[ — t, t] . Its top line is the object of interest, since it coincides with PNG. 




1 if there is a height line passing through (j, t) , 
if there is no such line , 



(18) 
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To write down the defining kernel for r\ T we set, as operators on l 2 = 



HMJ) = -ip(j-l)-ip(j + l), (19) 

hMj) = -V0'-i)-V'(j + i) + -V'(j) ) (20) 

T 

and B T the spectral projection on {H T < 0}. B T is known as discrete Bessel 
kernel. Then 

R T (j,t-j',t')=(e- tH *(B T -m(t-t')y H A . (21) 

The moments of f] T (j, t) are given by the formula analogous to (jSJ)- For large 
r, h (x,r) = 2y/r 2 — x 2 , \x\ < r. Thus i] T (x,t) is not stationary in t, which 
is reflected in (121}) by the fact that ^ H T . 

The correct edge scaling can be guessed as for Dyson's Brownian motion, 
where in spirit r is equated with N. We focus our attention on a space-time 
window of width r 2//3 and height r 1//3 centred at t = and j = 2r. Properly 
rescaled, compare with (|12|). H T becomes 

H T ijj{x) = t 2/3 (—iP(x-t- 1/3 )-iIj(x + t- 1/3 ) + 2iIj(x) + t 1/3 (x/t)iIj(x)) , (22) 

which converges to the Airy operator H as r — ► oo. This argument overlooks 
that even in rescaled coordinates the lines have still a systematic curvature. 
The correct limit is thus, [■] denoting integer part, 

lim T 1/3 r] T ([2T + t 1/3 (x - t% r 2/3 t) = <f>(x, t) (23) 

with the Airy random field. In particular, the top line converges to the 
Airy process A(t). 

We summarise our discussion as 

Theorem 3.1. Let h(x, r) be the PNG model in the droplet geometry. Then, 
in the sense of weak convergence of finite- dimensional distributions, 

lim r~ 1/3 (h(tr 2/3 , r) - 2r) = A(t) - t 2 . (24) 

T— »00 

Proof. For t = this is the celebrated result of Baik, Deift, Johansson^] on 
the length of the longest increasing subsequence of a random permutation. 
They use orthogonal polynomials and Riemann-Hilbert techniques in their 
asymptotic analysis. Johansson |1 6 j develops an approach through Fredholm 
determinants, which is also the basis of the proof by Prahofer, Spohn[lO] for 
the time-extended case. □ 
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Theorem 13. II is stated for the reference point x = 0. An analogous limit 
holds for any other reference point x = ar, \a\ < 1. 

For a space-time discrete version of the PNG droplet Johansson |17j proves 
the analogue of Theorem 13. II including tightness. While the details are long, 
it is likely that Theorem 13.11 can be strengthed to weak convergence of path 
measures. As established by Johansson^H], for the Aztec diamond the border 
between the frozen and disordered zones is also governed by the Airy process. 

4 Extensions 

What is so special about the droplet? From the point of view of growth 
processes other initial conditions look more natural, e.g. the initially flat 
surface h(x, 0) = 0. The RSK dynamics can be implemented for any choice of 
initial data and set of nucleation events. However the determinantal property 
of the resulting line ensemble is fragile. 

We illustrate our point by three examples. 

(i) half droplet. In a discrete space-time setting this problem is studied re- 
cently by Sasamoto and Imamura[19j. The rule is to simply suppress all 
nucleation events for x < 0. In addition there is a source at x = with 
rate 7. For < 7 < 1 the additional mass is incorporated in the bulk with 
no change in the macroscopic profile. The height fluctuations at x = sat- 
isfy (JTJ) with F 2 replaced by F 4 , i.e. the distribution function of the largest 
eigenvalue of a GSE random matrix in the limit of large N • 7 = 1 is crit- 
ical and F 2 in ((TJ) is replaced by F\ [20], i.e. the distribution function of the 
largest eigenvalue of a GOE random matrix|2I]. For 7 > 1 the macroscopic 
profile develops a linear portion, starting at x — 0, at height (7 + 7 -1 )t and 
joining tangentially the profile for 7 = 0. The height fluctuations are then 
of order yfr and Gaussian. The RSK dynamics results in a line ensemble, 
which for fixed r has the following weight. At x = r the boundary condition 
is hj{r,r) = j. The lines have jump size one and are not allowed to cross. 
Under these constraints the jump points are uniformly Lebesgue distributed. 
If this construction would be extended to x — — r with the boundary condi- 
tion hj(—r, r) = j, the resulting line ensemble is the one of the PNG droplet. 
For the half-droplet the lines end at x = and there is an extra weight from 
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the source at the origin, which reads 

o 

exp[(log 7 ) Yl (MM - ^-i(0,r) - 1)] . (25) 

j=-oo 

For 7 = 1 the extra weight (|23|) equals 1 and the weight for the line ensem- 
ble can be written as a product of determinants. However, the correlation 
functions are not determinantal because of the free boundary at x — 0. By 
an identity of de Bruijn, gap probabilities are square roots of determinants, 
thus yielding GOE edge statistics at x = in the scaling limit. For 7 = the 
lines at x = are constrained as /i2j(0, r) — h 2 j-i(0, r) = 1, which resembles 
Kramers degeneracy, thus GSE. An asymptotic analysis is available only for 
7 = 0, 1 with the result that r 2 / 3 close to the origin the height statistics 
is governed by the largest eigenvalue of the transition ensemble governing the 
crossover from GOE, resp. GSE, at x = to GUE in the bulk. For 7 > 1, at 
x = the top line separates the distance (7 + 7 _1 — 2)r from /i_i(0, r) which 
remains located roughly at 2r, as for 7 < 1. This explains the Gaussian 
fluctuations for 7 > 1. 

(ii) stationary PNG. If initially the upward steps are Poisson distributed 
with density p+ and the downward steps with density p_, then PNG on the 
whole real line is a space-time stationary growth process, provided p+p- = 1. 
To determine h(0, r) it suffices to know the nucleation events in the back- 
ward light cone {(x, t)\0 < t < r, |x| < r — t}. In fact it also suffices 
to only know the nucleation events in the forward light cone of the origin, 
{(#,£) |0 < t < r, \x\ < r}. The reason is that for stationary PNG model the 
intersection points of the diagonals {x = ±i} with the world lines for steps 
are again Poisson distributed. Thus we arrive at a set-up rather similar to the 
PNG droplet. The additional feature is that there are sources of nucleation 
events at the boundaries ±r. The sources are Poisson in time with rates a + , 
resp. a_, a + > 0, a_ > 0, stationarity being equivalent to = 1. As be- 

fore, a line ensemble is generated through the RSK dynamics. At time r, the 
boundary conditions hj(±r, r) = j, j = —1, —2, hold. However, h (±r, r) 
can now take arbitrary positive integer values. They have geometric weight, 
i.e. the weight 

exp[(loga ± )/i (±T,r)] • (26) 

The step points are Lebesgue distributed, constrained by non-crossing and 
the boundary conditions at ±r. Let Q be the collection of all such lines. 
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By running RSK backwards in time, one will end up with hj(0, r) = j and 
h (0,r) = n for r > but sufficiently small. Q decomposes accordingly 
as Q = U n >ofl n . Q n has the total weight (a + ct-) n Z , Z Q = exp[t(a: + — 
a! 1 ) + t(a_ — QtZ )]. Stationary growth corresponds to the sector O carrying 
weight Wo and probability Po = Wq/Zq for the lines. However, only the line 
ensemble with weight W on Q is determinant al. Fortunately, one has the 
simple relationship 

M{h (0,r) = k}) = Z 1 W ({h (0,r) = k}) (27) 
= Z, 1 (W({/i (0, r) = k}) - a + a_W({/i (0, r) = k - 1})) . 

Stationary PNG is described by the spatial derivative of a determinantal 
process. (|27)l serves as a starting point for an asymptotic analysis, see Baik, 
Rains [22] and Prahofer, Spohn[23j. 

(hi) flat initial conditions. One sets h(x, 0) = and allows for nucleation 
events on the entire real line. For the RSK construction it is necessary to 
first restrict to the periodic box [—£, £]. Then hj(—£, r) = hj(£, r). As before, 
the lines are constrained on non-crossing and unit jumps. In addition they 
have to satisfy 

min (hj(x, r) — h-Ux, r)) = 1 . (28) 
\x\<e v ' 

An asymptotic analysis of the line ensemble with these constraints is not 
available. By completely different methods Baik, Rains [2*U] prove that, first 
taking £ — * oo, the height fluctuations at x = are distributed as the largest 
GOE eigenvalue, i.e. F 1 replaces F 2 in (JTJ). On this basis, a possible guess for 
the full height statistics, in the limit r —> oo, is Dyson's Brownian motion at 
(5 = 1. We will have to wait for the next ICMP congress for a confirmation 
(or a counter argument). 
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